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Abstract (800 characters maximum including - 150 words) 

We propose to use high-sensitivity BATSE hard X-ray observations to test the thick-target and electric field 
acceleration models of solar flares. We will compare the predictions made by these models with hard X-ray spectral 
observations obtained with BATSE and simultaneous soft X-ray Ca XDC emission observed with the Yohkoh Bragg 
Crystal Spectrometer (BCS). The increased sensitivities of the BATSE and BCS (relative to previous detectors) 
permits a renewed study of the relationship between heating and dynamical motions during the crucial rise phase of 
flares. With these observations, we will: 

1. investigate the ability of the thick-target model to explain the temporal evolution of hard X-ray emission relative 
to the soft X-ray blueshift during the earliest stages of the impulsive phase. 

2. search for evidence of electric-field acceleration as implied by temporal correlations between hard X-ray spectral 
breaks and the Ca XEX blueshift. 

The proposed study will utilise hard X-ray lightcurve and spectral measurements in the 10-100 keV energy range 
obtained with the BATSE Large Area Detectors (LAD). The DISCLA and CONT data will be the primary data 
products used in this analysis. 
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1. INTRODUCTION 


This report discusses results of work performed under the Phase 3 Compton Gamma-Ray Obser- 
vatory (CGRO) Guest Investigator Program. The goal of this work is to study solar flare models 
by comparing their predictions with simultaneous hard and soft X-ray observations. We have 
used hard X-ray observations from the CGRO Burst and Transient Source Experiment (BATSE) 
and soft X-ray observations from the Bragg Crystal Spectrometer (BCS) and Soft X-ray telescope 
(SXT) on the Japanese Yohkoh spacecraft. 

We investigated a model in which solar flares are heated by large-scale electric currents and their 
associated DC-electric fields. DC-electric fields have been shown to provide a viable mechanism for 
heating plasma and for accelerating particles in solar flares (Holman 1985, Winglee et. al. 1991). 
In particular, quasi-static DC-electric fields parallel to the coronal magnetic field can accelerate 
thermal electrons until a steady-state current is established. Since the collisional drag on the 
electrons decreases with increasing velocity, the electrons with velocities above a critical velocity 
will undergo runaway acceleration to super-thermal energies (Dreicer 1959). These runaway 
electrons can produce significant nonthermal hard X-ray emission via thick-target interactions 
(Holman, Kundu, and Kane 1989; Benka and Holman 1992). Electrons with velocities below the 
threshold for runaway acceleration continue to provide energy to the ambient plasma (via Joule 
heating) producing thermal soft X-ray emission. 

2. WORK PERFORMED AND RESULTS 

We tested the DC-electric field model by comparing its predictions with the nonthermal electron 
heating rate deduced from BATSE hard X-ray observations and the thermal heating rate derived 
from Yohkoh soft X-ray line emissions. To perform this test, we formulated the DC-electric field 
model in terms of observable parameters. We expressed the ratio of the electron runaway rate 
N (electrons s _1 ) to the thermal heating rate Q (ergs s _1 ) as a function of electron temperature 
T, electron density n, and electric field strength E: 

N/Q = 2.5 x 10 8 T 7 ~ 3/2 (E d /E) 19/8 x exp[-2 1 / 2 (E D /E) 1/2 - (1/4)(E D /E)J (1). 

where Ep ~ 7 x lO^ngTV -1 (volts cm -1 ) is the Dreicer electric field, wdiich is the electric field 
strength for which all thermal electrons undergo runaway acceleration. The numerical subscripts 
denote the order of each parameter (e.g., 7V = Tj 10 7 , etc). This ratio has the useful property 
that it is independent of the geometry of the current system (such as the thickness and number of 
current, channels). Given that T and n can be constrained by soft X-ray observations, equation (1) 
provides a useful diagnostic for determining E. 
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We applied this diagnostic by using BATSE LAD Continuous (CONT) data to infer N and 
BCS Ca XIX spectra to infer Q. The CONT data consist of 16 channels in the ~ 20 to > 400 keV 
range acquired at 2.048 s time resolution. This resolution is comparable to the 3 s time resolution 
of the BCS detector. We estimated N based on fits of nonthermal power-law spectra of the form 
I = a if - "*, where e is the photon energy (keV), a i is the amplitude, and 7 is the spectral index. 
For a power-law differential electron spectrum above a cutoff energy E c , these fitted parameters 
are related to the electron number flux by the formula (cf. Lin and Hudson 1976): 

Nobs ^ 3 x 10 33 ai (7 - 1 ) 2 £( 7 - 1/2, 1/2)E C ~ 1 (2), 

where B(.r, y) is the beta function. For the cutoff energy, we used the critical energy that sep- 
arates thermal and nonthermal electrons (Benka 1991. Benka and Holman 1992). Specifically, 
E c = m e v c 2 / 2 where v c is the threshold velocity above which thermal electrons exceed the fric- 
tional force and runaway. This threshold velocity is given by v c = (Ed/E) ' v e where v e is the 
electron thermal velocity. Figure 1 demonstrates a power-law fit applied to BATSE CONT spectra 
for a flare observed on 199 September 6. 

We estimated Q by equating the heating rate produced by current-carrying electrons with the 
heating rate required to sustain the observed soft X-ray emission. Thus, Q 0 bs — dU/dt + P r T P c , 
where U = 2>nkTV is the total thermal energy, and P r and P c are the conductive and radiative 
energy loss rates. The electron temperature and density were deduced from spectral fits to the 
BCS Ca XIX spectra. These fits effectively yield the emission measure EM — f n 2 dV which is 
related to the density via the loop geometry parameters (length L and cross-sectional area A c ) 
derived from Yohkoh SXT images. 

Substitution of the observed ratio (N/Q) 0 bs into the left hand side of equation (1), reduces it to 
a transcendental equation with the electric field strength E as the single unknown. This equation 
can be solved by least-squares or simple iterative techniques. We have solved for E in two flares 
that were observed jointly by CGRO and Yohkoh on 17 February 1992 and 9 September 1992. 
The electric field strength was found to vary from below 10 -5 volts cm -1 during the preflare phase 
to a maximum of 10 -4 volts cm -1 at the peak of the hard X-ray impulsive phase. The maximum 
field strength was sufficient to account for the number flux of nonthermal electrons required to 
produce the observed hard X-rays and - when applied over a loop scale length - sufficient to 
explain the maximum energies (~ 100 keV ) of nonthermal electrons. 

We have found that the variation in hard X-rays is strongly tied to the behavior of the ratio 
e = E/Ed ■ For example, a factor 5 increase in e from preflare to flare maximum corresponds to 
over a factor 100 increase in hard X-ray emission. This strong sensitivity is largely due to the 
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exponential dependence of the runaway rate on c (cf. equation 1). This sensitivity suggests that 
the relative behaviors of hard and soft X-rays can be physically understood in terms of variations 
in the electric field strength and/or Dreicer field. 

4. SUMMARY 

We have combined BATSE hard X-ray and Yohkoh soft X-ray observations to test a model in 
which electric currents and their associated electric fields are responsible for the respective heating 
and particle acceleration that occur in solar flares. We have found that the DC-electric field model 
can be used to explain the number of flux of accelerated nonthermal electrons implied by hard X- 
ray observations and the thermal energetics implied by soft X-ray observations. We have used the 
model to compute the strength and evolution of the coronal electric field and to derive important 
flare parameters such as the current density and the degree of current filamentation in the corona. 

Our calculations have confirmed similar results obtained from previous applications of the DC- 
electric field model. We have found that the density of the flare plasma is a critical parameter 
that controls the evolution of impulsive hard X-rays. In particular, a small increase in density can 
effectively quench the runaway acceleration process by increasing the frictional drag on thermal 
electrons. Similar effects of increasing density on the hard X-ray emission were also found by 
Tsuneta (1985) in Hinotori data. 

Further details and results of our analysis are published in Zarro, Mariska, and Dennis (1995) 
and Love et al. (1995). We are presently using the DC-electric field model to analyze additional 
flares that were observed simultaneously by CGRO and Yohkoh. 
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Energy (keV) 

Interval 12 

09:03:52.45- 09:04:00.64 

f_vth_bpow parameters: 1.000e-05, 1.000, 0.5843, 4.444, 400.0, 1.930 


Figure 1. The BATSE Large Area Detector Continuous (LAD/CONT) hard X-ray lightcurves (upper panel) for a 
flare observed on 1992 September 6. A single power-law fit (lower panel) to the CON'T spectra in the 25-120 keV 
range during one of the hard X-ray burst maxima. 
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TESTING THE DC-ELECTRIC FIELD MODEL IN A SOLAR FLARE OBSERVED 
BY YOHKOH AND THE COMPTON GAMMA-RAY OBSERVATORY 

D. M. Zarro, 1 J. T. Mariska, 2 and B. R. Dennis 3 * 

Received 1994 February 22; accepted 1994 September 2 

ABSTRACT 

We apply a DC-electric field model to the analysis of soft and hard X-ray observations of a solar flare 
observed by Yohkoh and the Compton Gamma-Ray Observatory CGRO on 1992 September 6. The flare was 
observed simultaneously in the soft X-ray Ca xix line by the Yokhoh Bragg Crystal Spectrometer (BCS) and 
in hard X-rays (>50 keV) by the CGRO Burst and Transient Spectrometer Experiment (BATSE). A strong 
stationary component of Ca xix emission was present at the start of impulsive hard X-ray emission indicating 
an extended phase of heating prior to the production of energetic nonthermal electrons. We interpret the pre- 
flare Ca xix emission as a signature of Joule heating by field-aligned currents. We relate the temporal varia- 
tion of impulsive hard X-ray emission to the rate of runaway electron acceleration by the DC-electric field 
associated with the current. We find that the initial rise in hard X-ray emission is consistent with electron 
acceleration by a DC-electric field that increased from a preflare value of < 10” 5 V cm -1 to ^(9 ± 1) x 10“ 5 
V cm -1 at the time of the first hard X-ray peak and then remained constant during the rest of the impulsive 
phase. We attribute the increase in electric field strength to the formation of a current sheet at the reconnec- 
tion point of two loop structures. The decrease in hard X-ray emission after flare maximum is consistent with 
a reduction in the number of runaway electrons due to an increase in coronal density produced by chromo- 
spheric evaporation. The increased density quenches the runaway process by enhancing collisional thermaliza- 
tion of electrons. To avoid the generation of an unrealistically large magnetic field, the flaring region must be 
highly filamented into >10 6 oppositely directed current channels of ~30 cm width with an initial preflare 
current of ^3 x 10 10 A per channel. 

Subject headings : MHD — Sun: flares — Sun: X-rays, gamma rays 


1. INTRODUCTION 

Electric currents and their associated electric fields have 
been shown to provide a viable mechanism for heating plasma 
and for accelerating particles in solar flares (Tsuneta 1985; 
Holman 1985; Moghaddam-Taaheri & Goertz 1990). In par- 
ticular, quasi-static DC-electric fields parallel to the coronal 
magnetic field can accelerate thermal electrons until a steady 
state current is established. Since the collisional drag on the 
electrons decreases with increasing velocity, the electrons with 
velocities above a critical velocity will undergo runaway accel- 
eration to superthermal energies (Dreicer 1959; Spicer 1982). 
These runaway electrons can produce significant nonthermal 
hard X-ray (HXR) emission via thick-target interactions 
(Holman 1985; Holman, Kundu, & Kane 1989). Electrons with 
velocities below the threshold for runaway acceleration con- 
tinue to heat the ambient plasma (via Joule heating) thereby 
producing thermal soft X-ray emission. Hence, an attractive 
feature of electric fields is the potential for simultaneously 
explaining thermal and nonthermal processes occurring in 
solar flares. 

The mechanism for producing coronal electric fields is a 
subject of considerable debate. It is widely believed that such 
fields are generated by magnetic instabilities and reconnection 
processes that are involved in the triggering of a solar flare (cf. 
Tandberg-Hanssen & Emslie 1988). Numerical studies suggest 
a scenario in which electric field structures may be produced in 

1 Applied Research Corp., Solar Data Analysis Center, Code 682 NASA 
GSFC, Greenbelt, MD 20771. 

2 Naval Research Laboratory, Code 7673 Washington, DC 20375. 

3 Laboratory for Astronomy and Solar Physics, Code 682.2 NASA/GSFC, 

Greenbelt, MD 20771. 


coronal magnetic field configurations that are sheared progres- 
sively by photospheric footpoint displacements (Roumeliotis, 
Sturrock, & Antiochos 1994). This physical connection 
between loop current systems and flare phenonema is rein- 
forced observationally by high-quality vector magnetograph 
data that have revealed a close spatial association between 
vertical current components and Ha electron precipitation 
sites in the chromosphere (Canfield et al. 1993). 

In this paper, we present new results obtained using a DC- 
electric field model to interpret simultaneous soft and hard 
X-ray observations of a flare obtained by instruments onboard 
the Japanese Yohkoh and Compton Gamma-Ray Observatory 
(CGRO) spacecraft. The increased spectral sensitivities of these 
instruments (relative to similar detectors on preceding space- 
craft such as SMM) permit a study of the crucial early rise 
phase of the flare impulsive phase when signatures of plasma 
heating and particle acceleration are expected to be most 
revealing. 

2. OBSERVATIONS 

A GOES-c lass M3. 3 flare occurred at 08:57 UT on 1992 
September 6 in Active Region AR 7270 located at Sll W38. 
The flare was one of a series of events in this active region that 
was well-observed by Yohkoh (Shibata et al. 1993). It was 
observed simultaneously in the soft X-ray Ca xix (3.17 A) line 
by the Yohkoh Bragg Crystal Spectrometer (BCS) and in hard 
X-rays by the CGRO Burst and Transient Source Experiment 
(BATSE). 

The BATSE instrument (Fishman et al. 1989) consists of 
eight detector modules. Each module contains two Nal scintil- 
lation detectors: a Large Area Detector (LAD) and a Spectros- 
copy Detector (SD). Figure 1 a shows the evolution of the 
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Fig. 1.— (a) The BATSE Discriminator Large Area (DISCLA) hard X-ray 
count rate for the 50-100 keV channel plotted at 1.024 s time resolution. The 
vertical dotted lines mark the accumulation interval of the BCS Ca XIX 
preflare spectrum plotted in Fig. 2. 


Discriminator Large Area (DISCLA) hard X-ray count rate 
(1.024 s temporal resolution) in the 50-100 keV channel of the 
most sunward-facing LAD. Significant HXR emission was first 
detected at 09:02 UT. The emission reached a maximum at 
09:03:31 UT and was followed by three successive bursts each 

15-20 sin duration. x „ c 

The BCS (Culhane et al. 1991) is an uncolhmated full-Sun 
field-of-view spectrometer that consists of four detector chan- 
nels that span narrow wavelength ranges centered on emission 


1 . 5 * 1 0 5 ~~ 
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lines of highly ionized Fe xxvi, Fe xxv, Ca xix, and S xv. The 
nominal temporal resolution of the spectra is 3 s. Figure lb 
shows the light curve of the total count rate in the Ca xix 
channel. Enhanced Ca xix emission was observed 2-3 minutes 
before the peak of hard X-rays, with a strong stationary reson- 
ance line component (Fig. 2) already present at 09:01:57 
UT-the approximate start time of hard X-rays. 

2.1. Flare Geometry 

We use 119 Jim Be filter (centered on as 10 A) partial frame 
images from the Yohkoh Soft X-ray Telescope (SXT) to infer 
the geometry of the flare region. Figure 3 (Plate 11) shows a 
sequence of the available Be images during the interval 
09:02:24-09:03:56 UT, corresponding to the rise phase of 
hard X-rays. To facilitate relative comparisons, the images 
have been renormalized to the peak intensity of the last image 
in the sequence. These images indicate that the soft X-ray emis- 
sion originated in an approximately symmetric east-west (E- 
W) oriented looplike structure that was brightest initially at its 
center. At ~09:03:04 UT, the western leg of the loop bright- 
ened and eventually dominated the total emission from the 
loop structure. Hard X-ray images ( > 50 keV) obtained with 
the Yohkoh Hard X-ray Telescope (HXT) indicate that the 
impulsive phase burst emission occurred at the footpoints 
(denoted A and B) of the E-W loop structure (Yaji et al. 1993). 

The spatial and temporal sequences of soft and hard X-ray 
brightenings suggest a geometry in which the smaller E-W soft 
X-ray loop interacted with a larger overlying loop structure 
(with footpoints denoted by C and B in Fig. 3). This situation is 
similar to the emerging flux model discussed by Heyvaerts, 
Priest, & Rust (1977). Accordingly, we identify the apex of the 
smaller loop as the energy release site at which electrons were 
initially accelerated and subsequently propagated to the foot- 
points to produce thick-target hard X-ray emission. The 
separation of the footpoints A and B is approximately d * 
2875. Assuming a semicircular geometry, we deduce a corre- 
sponding half-length of L obs = 7td/2 ~ 3.2 x 10 cm for the 
E-W loop. We estimate the cross-sectional area of the E-W 
loop by computing the average width w along the region of soft 
X-ray emission connecting the A and B footpoints. Assuming a 
cylindrical loop volume and equating w with the loop diam- 
eter, we infer a cross-sectional area A obs = nw 2 /4 as 2.0 x 10 n 
cm 2 and a total volume V obi = 2A obi L obs a: 13 x 10 2 cm 



o 
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Fig. 2.— The BCS Ca XIX spectrum (15 s accumulation) at 09:01 :57 UT. 
the approximate time when the first significant increase in BATSE hard X-ray 
emission above 50 keV was detected. The smooth solid line superposed on the 
spectrum is a two-component synthetic fit to the Ca xix resonance and satellite 
lines The dotted line indicates the stationary component. The underlying 
smooth line is the blueshifted component. The wavelength scale of the Ca xix 
spectrum was calibrated by centering the resonance line during the late decay 
phase of the flare at the laboratory rest wavelength of 3.177 A 


2.2. Flare Energetics 

We determine the electron temperature T e of the stationary 
soft X-ray plasma from the relative intensities of the Ca xix 
resonance and dielectronic satellite lines (Bely-Dubau et al. 
1982). The line strengths are computed using a least-squares 
technique to fit synthetic profiles to the Ca xix spectra (cf. 
Fludra et al. 1989). Figure 2 shows a fitted synthetic spectrum 
superposed on the observed preflare spectrum at 09:01 : 57 UT 
The spectral fit implies a preflare temperature T e = (13 ± 2) 
x 10 6 K and a corresponding preflare emission measure of 
EM = j n 2 dV = (5 ± 2) x 10 47 cm' 3 . Using the loop volume 
estimated Vrom the SXT images and assuming a unity filling 
factor (see § 3), we derive the characteristic density n, a. 
[(EM F ob ,)]' 1 of the loop. The evolutions of T e and n e during 
the impulsive phase are plotted in Figures 4 a and 4b, respec- 

The Ca xix spectra for this flare were exceptional because of 
a strong blue asymmetry that was evident before the peak of 



890 


ZARRO, MARISKA, & DENNIS 


Vol. 44# 





Fig. 4. — (a) The electron temperature of the soft X-ray plasma deduced 
from fits of synthetic spectra to the BCS Ca xix resonance and satellite line 
profiles. To improve the signal- to- noise of each fit, the observed spectra were 
accumulated into 9 s intervals, (h) The electron density of the loop soft X-ray 
plasma deduced from the Ca xix emission measure assuming a simple loop 
geometry and a unity filling factor, (c) The iine-of-sight blueshift velocity 
deduced from the separation between the primary Ca xix resonance line and 
the secondary component. Note the first velocity peak at 09:02 UT (before 
significant HXR emission) and the secondary peak at 09:03: 10 UT (during the 
HXR rise phase). 

hard X-rays. We have modeled the asymmetry by including a 
blueshifted component in the above least-squares fitting 
analysis. Figure 2 shows the resulting two-component fit to the 
preflare spectrum at 09:01:57 UT. The separation of the sta- 
tionary and blueshifted components indicates an initial line-of- 
sight upflow velocity of v % 270 km s' 1 at the time of the 
preflare spectrum. Figure 4c shows the evolution of the blue- 
shift velocity during the impulsive phase. The upflow shows a 
secondary peak at 09:03 : 10 UT. To within the 9 s resolution of 
the BCS, this peak corresponds to the time when the SXT E-W 
loop feature showed enhanced brightening in the western leg. 
This temporal association suggests that the loop brightening 
was produced by gradual filling of the loop leg by chromo- 
spheric evaporation of soft X-ray plasma from the footpoint 
closest to the interaction site. 

We compute the thermal heating rate Q obs needed to sustain 
the soft X-ray emission from the loop by using the formula 

Gobs = dU/dt + P rad + P cond ergs s" 1 , (1) 

where U = 3 n e kT e V oh% is the total thermal energy of the 
plasma, P cond ^ 10“ 6 Tj /2 (2 A obs/ L ob J is the classical conduc- 
tive loss rate, and P rad * 1.5 x 10" l9 n 2 e T~ 1 2 K obs is the radi- 
ative loss rate (Antonucci, Gabriel, & Dennis 1984). The factor 
of 2 accounts for the assumed symmetry of the loop. Figure 5a 
shows the temporal variations of each of the terms in equation 
(1). Figure 5b shows the corresponding variation of Q obs . Note 
that the radiative loss term is negligible relative to the conduc- 
tive term during the flare. For sufficiently strong heating, the 
resulting high-temperature gradients can produce free- 


streaming electrons that carry the conductive flux. In this case, 
the conductive flux is given by the saturated flux F M * 
i n € m e v*, where m e is the electron mass and v e = (kT/m e ) i/2 is 
the electron thermal velocity (Manheimer & Klein 1975). The 
variation of P Mt = 2^ ob5 F Mt is shown in Figure 5a. Since the 
classical conductive term is below the saturated value for the 
adopted values of T e and n € , we use the classical term in equa- 
tion (1). An additional limitation on the heat flux arises when 
bulk streaming of the electrons sets up ion-acoustic waves 
which enhance the collision frequency. In this case, the conduc- 
tive flux is given by the anomalously limited heat flux 
(Tandberg-Hanssen & Emslie 1988). In § 3, we find that this 
anomalous-heating effect is not expected to be significant for 
the flare under study. 


3. ANALYSIS 

The 1992 September 6 flare is typical of many flares 
observed by Y ohkoh in that it showed evidence of strong high- 
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Fig. 5 -Ha) The computed time rate of change of the thermal energy dV/dt 
in a loop with uniform temperature and density derived from BCS Ca xix 
observations. The classical conductive P M and radiative P ^ energy-loss 
rates, as well as the saturated conductive loss rate P m are shown for compari- 
son. (b) The total heating rate of the soft X-ray plasma Q ^ = dVjdt + f> 

+ ?tu i- For the Ca xix-derived temperatures and densities, the heating rate Is 
dominated by the classical conductive loss term. 
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temperature Ca xix emission before the onset of impulsive 
hard X-rays (Mariska et al. 1993). Again, as is typical in many 
flares, this stationary component was accompanied by a 
weaker blueshifted component which indicates evaporative 
upflows that commenced prior to the HXR impulsive phase. 
Similar behaviors were found in flares observed by SMM 
(Doschek et al. 1986; McClements & Alexander 1989). These 
observations place severe constraints on models of solar flares. 
For example, the thick-target driven chromospheric evapo- 
ration model predicts that the blueshifted component will ini- 
tially dominate the Ca xix profile when the plasma is being 
heated and evaporated by energetic electrons (Mariska, 
Emslie, & Li 1989; Li, Emslie, & Mariska 1990). 

The presence of a dominant preflare stationary Ca xix com- 
ponent appears to contradict the thick-target electron beam 
model unless some additional form of preflare heating is 
invoked (Mariska & Zarro 1991). If this heating is purely 
thermal in origin, then equation (1) implies a plasma heating 
rate of Q obt = 2 x 10 28 ergs s' 1 before any significant hard 
X-ray emission was observed in the 1992 September 6 flare. 
The cause of this preflare heating is controversial. The various 
explanations that have been postulated include the following: 
proton beams that heat the chromosphere without producing a 
direct HXR signature (Simnett & Haines 1990); anomalous 
conductive heating by thermal electrons (Takakura 1992); and 
sudden pinch-heating at the loop top (Feldman et al. 1994). 

We consider an alternative explanation whereby the preflare 
soft X-ray heating is produced by Ohmic dissipation in thin 
current channels (or sheets) that are aligned parallel to the 
direction of the magnetic field in the E-W loop. These current 
channels generate an associated DC-electric field that simulta- 
neously accelerates electrons along the magnetic field to 
produce impulsive hard X-rays. In the following analysis, we 
will combine the soft and hard X-ray observations of the 1992 
September 6 flare to derive how the DC-electric field strength 
must have varied in order for this explanation to be plausible. 

Holman, Kundu, & Kane (1989) give the following formulae 
for the total Joule heating rate Q and electron acceleration rate 
S in a field-aligned current system with magneticjeld strength 
B G and associated electric field strength £ V cm ~ 1 

Q= 1.11 x 10 2l n,.4 18 B 2 T\ l2 v e (E/E D ) ergs s~ 1 , (2a) 

N ~ 2.83 x 10 2 V^i8 B 2 7'7 1/2 v ,(£ d /£) 11/8 

x exp [ — (2£ D /£) W2 - (£ D /4£)> ■"* , (2b) 

where v, = 10n 9 T 7 3/2 s~ l is the thermal collision frequency, 
n f is the number of current channels, 4 lg cm 2 is the area of 
individual channels, and E D = 7 x 10~ 6 n 9 T ? 1 V cm 1 is the 
Dreicer electric field. The numerical subscripts denote the 
order of each parameter (e.g., T 7 = TJ10 7 , etc). For £ < £d> 
the only electrons that runaway are those with velocities above 
a critical threshold velocity v c = {E 0 jE) ul v,. In particular, for 
£/£„ < 0.01, Joule heating will dominate over acceleration 
since there are few electrons with velocities above this critical 
threshold. For increasing values of £/£ D , acceleration will 
dominate as a larger fraction of the electrons enter the 
runaway regime and produce increasing nonthermal HXR 
emission. Ultimately, all electrons will runaway when £ = E D . 
This dependence of increasing HXR emission on increasing 
runaway rate suggests an empirical relationship F HX r = <*N , 
where F HX r is the intensity of nonthermal HXR emission pro- 
duced by the runaway electrons and a is a constant of pro- 
portionality that is to determined from the observations. 


The relative importance of acceleration versus heating can 
be parameterized by the ratio N/Q, which has the useful pro- 
perty that it depends only on T e , n e , and £. In particular, given 
the Ca xix-derived measurements of T e and n e , this ratio pro- 
vides a novel diagnostic for the unknown electric field strength. 
Accordingly, we can write 

F HX JQ ob , = xN/Q , (3) 

where the right-hand term is the theoretical ratio from equa- 
tions (2a) and (2b), and the left-hand term is the ratio of the 
observed nonthermal HXR emission to the thermal heating 
rate implied by BCS soft X-ray observations. We assume that a 
is constant during the flare and set F HXR equal to the total 
BATSE count rate intensity summed above 50 keV with the 
preflare background subtracted. (The emission below 50 keV is 
excluded in order to minimize possible contamination by 
thermal emission.) For the temporal variation of £, we assume 
a simple functional form in which E increased during the HXR 
rise phase and remained constant after the first peak of hard 
X-rays. Using T e and n e from the soft X-ray observations, we 
solve first for the constant values of E and oc after this first peak. 
We use a two-parameter least-squares fit of equation (3) to 
F H xr during the impulsive phase of the flare when BATSE 
showed the strongest burst emission (i.e., from the time of the 
first HXR peak to 09:06 UT). Having solved for the constant a 
after the peak, we perform a backward solution of equation (3) 
to derive the temporal variation in E from preflare to HXR 
maximum. 

Figure 6 a shows the results of the fitting process. The com- 
puted F hxr is superposed on the BATSE intensity above 50 
keV. The computed £ during the preflare and impulsive phase 
is shown in Figure 6b. We find that the overall evolution of 
hard X-rays during the impulsive phase is consistent with a 
DC-electric field strength that increased montonically from 
a preflare value of <10“ 5 V cm -1 to a value of 
£ _ (9 + i) x 10“ 5 V cm -1 at the time of the first peak of the 
HXR burst and remained constant thereafter. The evolution of 
the hard X-ray emission during the impulsive phase is tied 
strongly to the variation of the ratio e = £/£d- Figure 6b 
shows that a factor of 5 increase in € from preflare to HXR 
maximum corresponds to over a factor of 100 increase in HXR 
emission. This strong sensitivity is largely due to the exponen- 
tial dependence of the runaway rate on e (cf. eq. [2b]) which 
stems from the Maxwellian energy distribution function of the 
thermal electrons. 

The assumption of a unity filling factor in the calculation of 
n e implies that the loop volume is filled uniformly with soft 
X-ray emitting plasma that is at the temperature derived from 
BCS Ca xix. It is likely that the flare volume is filamented into 
finer structures with intrinsically different temperatures and 
densities. Since the full-Sun BCS instrument averages over the 
entire flaring active region, the Ca xix-derived values of T e and 
n e provide lower bounds to the actual temperature and density 
of finer structures within the loop. Given the computed value 
of £, we can estimate the filamentation of the DC-electric field 
current system and its relationship to the coronal filling factor. 
We use the time of the first HXR peak as the reference point of 
this calculation. Assuming a laminar sheet geometry, the total 
volume of n s current channels is given by V s = n s L s w s dr , where 
L s is the length, is the width, and Sr ( « u s ) is the thickness of 
individual sheets. We derive Sr from the constraint that the 
self-induction magnetic field of individual current channels 
must be less than the ambient coronal magnetic field strength. 
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Holman (1985; eq. [8a]) gives 

Sr £ 808 B 2 rtg 1 Tf xn E 0 !E cm , (4) 

where B 2 is the coronal magnetic field strength (in 100 G units) 
Substituting the values of n g , T ? , and E D /E at the first HXR 
peak gives Sr $ 30 B 2 cm and V s < 30n, L, w, B 2 . The assump- 
tion of a unity filling factor implies that the total assembly of 
sheets uniformly fills the loop volume such that V % — V 
Adopting a typical coronal field value of 100 G and assuming 
that the channels extend along the full length of the loop (with 
w J and L, equal to the total loop length 2 L ob< = 6 x 10 9 cm), 
^derive n, > 4 x 10 6 and <5r < 30 cm. The observation of 
HXR emission at opposite footpoints of the coronal loop sug- 
gests that a fraction (~50%) of the sheets have field-aligned 
currents that are directed in opposite directions along the loop. 

The ability to observe such thin structures is, of course, 
beyond the resolution capability of current or proposed instru- 
mentation. However, if we assume that the same number of 
individual sheets existed before the impulsive phase, then we 
can estimate the preflare current /, within each sheet that is 
required to explain the observed soft X-ray emission by Joule 
heating. We use the formula Q ob , = J£F obI , where the right- 
hand term is the total Joule heating rate and J is the current 
density. Substituting the values of £ and Q obt corresponding to 
the preflare time 09:01:57 UT, we derive J = 4.5 x 10 8 
statamp cm ~ 2 or, equivalently, J = 0.15 A cm 2 . Hence, the 
current in a single sheet with cross-sectional area w,Sr is given 


by l, Jw j dr ~ 2.7 x 10*° A. Invoking a smaller fling factor 
( < 1) would imply a larger density for the soft X-ay plasm 
within the sheets and, hence, a greater frictional toag oo the 
thermal electrons. Consequently, stronger electric fclds (ani 
hence, stronger current heating) would be required m older to 
produce the same nonthermal HXR emission by raaaway elec- 
trons. An additional consequence of increased (falsity is to 
elevate the radiative cooling P r>d . This enhanced radiatm 
cooling is crucial since it can potentially counterbalance the 
increased Joule heating by the coronal current, tfas restoring 
the overall heating rate to remain consistent with toe observed 
thermal soft X-ray emission levels. 

Finally, we compare our results with those deriwl by Benka 
(1991) and Benka & Holman (1992). The latter invotigators fit 
hybrid thermal/nonthermal electron distribution fcnetions to 
the high-resolution germanium HXR spectra of Lin et aL 
(1981). Thermal fits to the low-energy component of these 
spectra ^indicate “super hot ” temperatures of (30-Hto) x 10 6 K 
in ~ 10 current sheets with typical maximum field strengths dt 
9 . 9 x 10 Vcm *. Such temperatures are considecably hotter 
than the peak Ca xix— inferred temperatures ( ~2D x 10 6 IQ 
used in our analysis. A direct consequence of ligher tem- 
peratures is that more electrons are available far runaway 
acceleration and, hence, weaker electric fields aad/or fewer 
current sheets are required to produce the same ■onthermrf 
HXR emission. 

4. SUMMARY AND CONCLUSIONS 

We have explored a new method for analyzing simultaneous 
soft and hard X-ray observations of solar flares. The method is 
based on the assumption that the temporal evolution of impul- 
sive hard X-ray emission is related physically to the runaway 
acceleration rate of electrons in a DC-electric field (Tsuneta 
1985; Holman 1985). The ratio of the electron runaway rate to 
the rate of Joule heating by a field-aligned current is a straight- 
forward function of the temperature, density and dectric fieW 
strength in the plasma (Holman, Kundu, & Kane 1989). Using 
Ca xix soft X-ray spectra to compute the temperature and 
density, we have solved for the strength and variation of the 
electric field that is required to simultaneously explain the 
observed preflare soft X-ray heating and the impulsive phase 
hard X-ray emission. 

We have applied this technique to the flare of 1992 Septem- 
ber 6 and developed the following scenario to explain the 
sequence of hard and soft X-ray observations: 

1* Prior to the impulsive phase, the ratio of the DC-electric 
field to the Dreicer field is sufficiently low that electron 
runaway acceleration is negligible compared to Joule heating 
by the coronal currents. Consequently, thermal soft X-ray 
emission dominates over nonthermal HXR emission. The 
Joule heating raises the temperature of the loop to produce the 
enhanced soft X-ray brightening at the loop top and the domi- 
nant preflare Ca xix stationary component. Thermal conduc- 
tion from the loop top transports energy to the footpoints 
where it is converted to mass motions, thus, accounting for the 
preflare soft X-ray Ca xix blueshifts. 

2. During the interval 09:02-09:03 UT, the DC-electric 
field strength increases. Based on high spatial resolution soft 
X-ray images, we postulate that this increase is triggered by 
magnetic reconnection between two intersecting loops. Recent 
work by Tsuneta (1993) indicates that magnetic reconnection 
can enhance the field-aligned coronal current and associated 
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DC-electric field by driving turbulent convective motions in 
the current sheet that is formed at the reconnection site. We 
note that the maximum value of 6 = 0.18 reached during this 
flare is below the level at which turbulence-induced instabilities 
(and associated anomalous heating effects) are expected to 
become pronounced (Tsuenta 1985). The increasing £ raises 
the Joule heating rate, thereby increasing the loop temperature 
and shifting a greater population of thermal electrons closer to 
the critical threshold for runaway acceleration. Consequently, 
an impulsive increase in nonthermal HXR emission is produc- 
ed as these electrons enter the runaway regime and are acceler- 
ated by the enhanced electric field. The bulk of the HXR 
emission is produced by thick-target interactions as the accel- 
erated electrons propagate along the loop magnetic field from 
the reconnection site to the chromospheric footpoints. The 
initial energy deposited by the electrons during the first HXR 
peak heats the footpoints and fuels additional chromospheric 
evaporation. This additional evaporation accounts possibly for 
the secondary peak in the Ca xix blueshift velocity at 09:03: 10 
UT (Fig. 4c). 

3. The density of the loop increases rapidly (Fig. 4b) during 
the impulsive phase as thick-target heating by the runaway 
electrons and thermal conduction from the Joule-heated loop 
plasma continue to drive chromospheric evaporation. The 
increasing density results in a higher value for the Dreicer field. 
Since £ is assumed constant after the first HXR peak, the 
increase in E D produces a decrease in e after this peak and a 
corresponding reduction in the rate of electron runaways. 
Physically, the effect of increasing n e is to amplify the frictional 
drag on thermal electrons thereby inhibiting runaway acceler- 
ation (Holman 1985). The net result is a quenching of the 
nonthermal HXR emission. The consequences of increasing n e 
on the HXR behavior were noted by Tsuneta (1985). Applying 
a DC-electric field model to a flare observed by Hinotori , he 
interpreted the overall evolution of HXR emission as a tran- 
sition from particle acceleration dominating under low-density 
conditions early in the impulsive phase to Joule heating 
dominating in the higher density conditions prevailing later in 
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the flare. He attributed the increase in density to chromo- 
spheric evaporation. The present Ca xix blueshift observations 
support this scenario. 

Our analysis has shown that the relative behaviors of hard 
and soft X-rays can be understood in terms of simultaneous 
heating and acceleration in multiple current channels. With 
several simplifying assumptions, our application of the DC- 
electric field model has provided a novel diagnostic for measur- 
ing the strength and temporal variation of the electric field 
within these channels, and a method for determining some of 
the physical parameters of these channels. Crucial questions 
relating to how such channels are created, whether they are 
stable, and how the current system is closed in the lower atmo- 
sphere are beyond the scope of the present work. Equally 
important are the issues of the temperature of the plasma 
within the sheets and the contribution to the overall energy 
balance due to heating by the runaway electrons themselves. 
The temperature values inferred from Ca xix probably under- 
estimate the actual sheet temperature. A more realistic deter- 
mination can be made in future using BCS observations of Fe 
xxv-xxvi line transitions. The effect of electron heating can be 
addressed by relating the electron runaway acceleration rate to 
the number flux of electrons precipitating into the chromo- 
sphere (cf. Lin & Hudson 1976). In particular, spectra! fits to 
high-resolution BATSE spectra can be used to express the 
observed nonthermal HXR flux in terms of the electron beam 
parameters, thus, eliminating the proportionality constant a 
and permitting a more direct solution for the electric field 
strength (Zarro et al. 1995). 
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Fig. 3. — A scries of Yohkoh SXT partial frame Bc-filter images at 10 s intervals spanning the initial rise phase of impulsive hard X-rays (09:02: 24-09 .03 : 56 UT). 
North is up and West is to the right. The images have been normalized to the same exposure times (1 s) and pixel sizes (2.45 arcsec). Each panel consists of 64 x 64 
pixels, corresponding to 157" x 157" square. To facilitate relative comparisons, the intensity of each pixel in each image has been renormalized to the intensity of the 
corresponding pixel in the last image of the sequence at 09:03:56 UT. The sequence of brightenings suggests the interaction of loops: a smaller east-west oriented 
loop (with footpoints marked A and B) connecting with a larger overlying loop (with footpoints marked C and B). The bright kernel of emission between A and B 
marks the interaction site. The western legs of the interacting loops show a gradual brightening that commenced at ^09:02:44 UT and increased monotonically 
during the flare, eventually dominating the total soft X-ray emission from the two loops. This brightening indicates chromospheric evaporation of plasma from the 
common footpoint B of the two loops. 
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